A measurement of lifetime differences in the 
neutral D-meson system 

The FOCUS Collaboration 

J.M. Link a , V.S. Paolone 8 " 1 , M. Reyes a ' 2 , P.M. Yager a , 
J.C. Anjos b , I. Bediaga b , C. G6bel b < 3 , J. Magnin b ' 4 , 
J.M. de Miranda b , I.M. Pepe b ' 5 , A.C. dos Reis b , 
F.R.A. Simao b , M.A. Vale b , S. Carrillo 6 , E. Casimiro 6 ' 6 , 
H. Mendez 6 ' 7 , A. Sanchez-Hernandez , C. Uribe 6 ' 8 , 

F. Vazquez , L. Cinquini d ' 9 , J. P. Cumalat d , J.E. Ramirez d , 

B. O'Reilly d , E.W. Vaandering d , J.N. Butler 6 , 
H.W.K. Cheung 6 , I. Gaines 6 , P.H. Garbincius 6 , L.A. Garren 6 , 
E. Gottschalk 6 , S.A. Gourlay 6 ' 10 , P.H. Kasper 6 , 
A.E. Kreymer 6 , R. Kutschke 6 , S. Bianco f , F.L. Fabbri f , 
S. Sarwar f , A. Zallo f , C. Cawmeld 8 , D.Y. Kim g , K.S. Park 8 , 
A. Rahimis, J. Wiss g , R. Gardner 11 , Y.S. Chung 1 , J.S. Kang 1 , 
B.R. Ko 1 , J.W. Kwak 1 , K.B. Lee 1 , S.S. Myung 1 , H. Park 1 , 

G. Alimontij, M. Boschim-i, D. Brambillaj, B. Caccianiga-', 
A. Calandrino j , P. D'Angelo j , M. DiCorato j , P. Dini j , 

M. Giammarchi J , P. Inzani J , F. Leveraro J , S. Malvezzi J , 
D. Menascej, M. Mezzadri-i, L. Milazzo^, L. Moroni 
D. PedriniJ, F. Prelz j , M. Rovere j , A. Sala j , S. Sala j , 
T.F. Davenport III k , V. Arena , G. Boca^, G. Bonomi^' 11 , 
G. Gianini^, G. Liguori^, M. Merlo^, D. Pantea^ 12 , S.P. Ratti £ , 
C. Riccardi^, P. Torre , L. Viola , P. Vituk/, H. Hernandez 111 , 
A.M. Lopez 111 , L. Mendez m , A. Mirles m , E. Montiel m , 
D. Olaya 111 ' 13 , J. Quinones 111 , C. Rivera 111 , Y. Zhang 111 ' 14 , 
N. Copty 11 ' 15 , M. Purohit 11 , J.R. Wilson 11 , K. Cho°, 
T. Handler , D. Engh p , W.E. Johns p , M. Hosack p , 
M.S. NehringP' 16 , M. Sales p , P.D. Sheldon p , K. Stenson p , 
M.S. Webster p , M. Sheaff q , Y.J. Kwon r 

^University of California, Davis, CA 95616 
h Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brazil 



Preprint submitted to Elsevier Preprint 



7 February 2008 



C CINVESTAV, 07000 Mexico City, DF, Mexico 
d University of Colorado, Boulder, CO 80309 
e Fermi National Accelerator Laboratory, Batavia, IL 60510 
i Laboratori Nazionali di Frascati delVINFN, Frascati, Italy, 1-00044 
^University of Illinois, Urbana- Champaign, IL 61801 
^Indiana University, Bloomington, IN 47405 
1 Korea University, Seoul, Korea 136-701 
} INFN and University of Milano, Milano, Italy 
^University of North Carolina, Asheville, NC 28804 
Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy 
m University of Puerto Rico, Mayaguez, PR 00681 
n University of South Carolina, Columbia, SC 29208 
° University of Tennessee, Knoxville, TN 37996 
p Vanderbilt University, Nashville, TN 37235 
q University of Wisconsin, Madison, WI 53706 
r Yonsei University, Seoul, Korea 120-749 



Abstract 

Using a high statistics sample of photoproduced charm particles from the FOCUS 
experiment at Fermilab, we compare the lifetimes of neutral D mesons decaying via 
D° — > K~tt + and K~ K + to measure the lifetime differences between CP even and 
CP odd final states. These measurements bear on the phenomenology of D° — D° 
mixing. If the D° — > K~ir + is an equal mixture of CP even and CP odd eigen- 
states, we measure y C p = (r(CP even)-r(CP odd))/(r(CP even)+r(CP odd)) = 
0.0342 ± 0.0139 ± 0.0074. 
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This paper contains a comparison of the lifetime of a CP even final state, 
D° — > K~K + , to the lifetime of a CP mixed decay, D° — > K~ir + . The lifetime 
measurements are made using high signal-to-background D° samples consist- 
ing of 10 331 decays into K~K + , and 119 738 decays into K~tt + . Throughout 
this paper, unless stated explicitly, the charge conjugate is implied when a 
decay mode of a specific charge is stated. 

If CP violation in neutral D-meson decays is negligible, the even CP and odd 
CP combinations of the D° and D° are mass eigenstates with well defined 
masses and widths. To the extent that D° <-> D° mixing transitions occur, 
the widths of the CP even and odd states may differ. This paper reports a 
new, direct measurement of yep = (r(CP even) — T(CP odd))/(r(CP even) + 
T(CP odd)). Throughout this paper, we refer to the width asymmetry between 
neutral D CP even and odd eigenstates as ycv to differentiate it from a mix- 
ing parameter (generally called y) which is the fractional width asymmetry 
between true mass eigenstates and could differ from j/cp to the extent that 
charm decays violate CP symmetry. 

Under the assumption that the decay D° — > K~n + is an equal CP even - 
odd mixture, the width difference asymmetry (yep) is related to the measured 
lifetimes via: 

r(CP even) - T(CP odd) _ t(D° -> K-n+) 
VCP ~ r(CP even) + T(CP odd) ~ t(D° K~K+) 

Because D° — > K~n + is assumed to be a mixed state, a sizeable width dif- 
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ference between the CP even and odd lifetimes could, in principle, create a 
deviation from a pure exponential time evolution. Given the present limits on 
y [1], this deviation is safely ignored given the scale of our statistical precision. 
Therefore, we will fit both lifetimes assuming a pure exponential decay. 

The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. FOCUS is a 
considerably upgraded version of a previous experiment, E687 [2]. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. We obtained a sample of 
over 1 million fully reconstructed charm particles in the three major decay 
modes: D° — > K~n + , K~7t + 7t~tt + and D + — > K~n + 'K + . We briefly discuss 
those aspects of the detector which are particularly relevant for this analysis. 

The FOCUS detector is a large aperture, fixed-target spectrometer with ex- 
cellent vertexing and particle identification. A photon beam is derived from 
the bremsstrahlung of secondary electrons and positrons with an pa 300 GeV 
endpoint energy produced from the 800 GeV/c Tevatron proton beam. The 
charged particles which emerge from the target are tracked by two systems 
of silicon microvertex detectors. The upstream system, consisting of 4 planes 
(two views in 2 stations), is interleaved with the experimental target, while 
the other system lies downstream of the target and consists of twelve planes of 
microstrips arranged in three views. These detectors provide high resolution 
separation of primary (production) and secondary (decay) vertices with an 
average proper time resolution of ~ 30 fs for 2-track vertices. The momentum 
of a charged particle is determined by measuring its deflections in two anal- 
ysis magnets of opposite polarity with five stations of multiwire proportional 
chambers. Three multicell threshold Cerenkov counters are used to discrimi- 
nate between electrons, pions, kaons, and protons. Our combination of a high 
quality, high resolution vertex detector, and an excellent Cerenkov system al- 
lows us to obtain clean and copious charm meson samples without the need 
of D* tagging although we employ both a tagged and untagged sample in this 
analysis. 

Throughout this analysis we have chosen cuts designed to minimize non-charm 
backgrounds as well as reflection backgrounds from misidentified charm de- 
cays. To minimize potential systematic error, we use cuts and selection tech- 
niques which create very little bias in the reduced proper time. The reduced 
proper time is a traditional lifetime variable used in fixed-target experiments 
which use the detachment between the primary and secondary vertex as their 
principal tool in reducing non-charm background. The reduced proper time 
is defined by t' — (£ — Nai) / {f3 r jc) where I is the distance between the pri- 
mary and secondary vertex, <j£ is the resolution on £, and N is the minimum 
"detachment" cut required to tag the charm particle through its lifetime. If 
absorption and acceptance corrections are small enough that they can be ne- 
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glected, and if <Ji is independent of £, one can show that the t' distribution 
for decaying charmed particles, in the absence of mixing effects, will follow an 
exponential distribution. These assumptions are very nearly true in FOCUS. 

With a few important differences, many of the basic cuts and analysis algo- 
rithms are described in reference [3]. We will summarize all important analysis 
issues here as well. Both states (D° — > K~n + and K~~K + ) were obtained using 
a data set which required a minimum detachment of the secondary vertex from 
the primary vertex of 2.5 and a high quality secondary vertex with a con- 
fidence level exceeding 1%. The primary vertex was found using a candidate 
driven vertex finder where a primary vertex was found by intersecting (nucle- 
ating) tracks about a "seed track" constructed using the secondary vertex and 
the reconstructed D momentum vector. The candidate driven algorithm finds 
the primary vertex with relatively high efficiency even at very low detachment. 
As we will demonstrate later in Figure 2, the use of candidate driven vertex 
finder meant that there was essentially no time dependent efficiency correction 
required to fit the reduced proper time to an exponential form. As a result, 
systematic uncertainties are greatly reduced since reliance on the Monte Carlo 
used to compute efficiency corrections is very minimal. 

We begin by describing additional vertexing and kinematic cuts that were 
used to reduce the background to D° — > K~ir + and K~K + . These cuts, along 
with more stringent detachment cuts and particle identification cuts, were 
used to obtain our final sample. To maximize our yields while maintaining 
good signal-to-noise, we allow an event to enter our sample through either a 
D* tagged path or an inclusive path with additional clean-up cuts. Both paths 
required that the primary vertex fell within the boundaries of our segmented 
target. 

The tagged path admitted any candidate consistent with the decay D* + — > 
D°n + by virtue of having a D* — D° mass difference within 3 MeV/c 2 of nom- 
inal. The very powerful tagging cut reduced the background to the extent that 
no additional kinematic cuts were required. The inclusive path required that 
anticipated proper time resolution for a given event satisfied a if (Pjc) < 60 fs, 
and required that the two tracks did not have grossly asymmetric momenta 
(| Pi - P 2 \/(Pi + Pi) < 0.70). We found from Monte Carlo studies, that 
background from partially reconstructed charm decays were often very highly 
asymmetric and effectively eliminated by the momentum asymmetry cut. We 
found that the proper time resolution cut was extremely effective at eliminat- 
ing long lived backgrounds. Rather than separately fitting samples from the 
two paths, we combined the candidates admitted through either path into a 
single sample (in a way that insured no double counting) for the purpose of 
fitting. In doing this, we are implicitly assuming that the D° -> K~K + and 
D° — > K~ K + have equal lifetimes and hence there is no need to discriminate 
between these two decays. 
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Because misidentified Cabibbo-allowed decays can be a significant background 
to the suppressed process D° K~K + , we have studied the charm particle 
lifetimes using a variety of Cerenkov cuts. The Cerenkov particle identifica- 
tion cuts used in FOCUS are based on likelihood ratios between the various 
stable particle identification hypotheses. These likelihoods are computed for 
a given track from the observed firing response (on or off) of all cells within 
the track's (/3 = 1) Cerenkov cone for each of our three, multicell thresh- 
old Cerenkov counters with pion thresholds of 4.5, 8.4, and 17.4 GeV/c. The 
probability that a given track will fire a given cell is computed using Pois- 
son statistics based on the predicted number of photoelectrons striking the 
cell's phototube under each particle identification hypothesis and an inten- 
sity dependent accidental firing rate determined for each of the 300 cells. The 
product of all firing probabilities for all cells within the three Cerenkov cones 
produces a x 2 -like variable called Wi = —2 log(likelihood) where i ranges over 
the electron, pion, kaon, and proton hypotheses. 

An important Cerenkov cut, which we will call kaonicity, is defined as AWk = 
W w — Wk- A kaonicity cut of AWk > 4 implies that the track we are assigning 
to the kaon has an observed Cerenkov pattern under the kaon hypothesis that 
is favored over that of the pion hypothesis by a factor of e 2 = 7.39. Over 
the momentum spectrum of typical tracks in FOCUS, pions have significantly 
different Cerenkov response than kaons, and only a small fraction (typically < 
15%) have AWk > 0. Thus, even a mild cut on the likelihood ratio favoring the 
kaon hypothesis on the kaon candidate track reduces backgrounds to Cabibbo- 
favored decays such as D° — > K~ir + , K~it + it~it + and D + — > K~tt + it + by a 
factor of ~ 10. 

As we will describe later, in order to assess systematic uncertainty, we varied 
the kaonicity and detachment cuts used to obtain the final sample prior to 
lifetime fitting. For each fit, a single detachment and kaonicity cut is used 
for both the K~ir + and K~ K + signal obtained through either the tagged or 
inclusive path. To further reduce backgrounds for the inclusive path of the 
signal, we required for tracks reconstructed as kaons, that the proton light 
pattern is not favored over the kaon hypothesis by more than Wk — Wp = 3. 
The inclusive path also required a pion consistency cut for the "pion" track 
in K~tt + decay that demanded that no particle hypothesis is favored over the 
pion hypothesis with a AW exceeding 2. No additional Cerenkov cuts beyond 
kaonicity were used for the tagged path. 

The mass distributions of Figure 1 illustrate several of the signals used in 
this analysis and specifically the use of these likelihood-based Cerenkov cuts. 
Figure 1(a) shows the D° — > K~it + signal obtained after requiring a kaonicity 
cut of AWk > 4 for the K~ . Figures 1(b) and 1(c) illustrate the use of tight 
kaonicity cuts to reduce the D° — > K~n + reflection background to D° — > 
K~K + . The reduction of the K~ir + reflection is evident as AW K , applied to 
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both tracks, is raised from 1 to 4. The D° — > K~K + signal yield is estimated 
using a Gaussian signal peak over a background consisting of a 5th-order 
polynomial to represent general backgrounds along with a K~n + reflection 
line shape taken directly from Monte Carlo but scaled by a fit parameter to 
best match the data. 




Fig. 1. (a) Signal for D° — > K ir + with a detachment cut of £/a > 5 and 
W n -W K > 4. The yield is 119 738 K'tt+ signal events. 

Signals for D° — > K~ K + with a detachment cut of tja > 5. The reflection 
in the background at higher masses is due to contamination from misidentified 
D° -» K--K+. (b) Requiring W V -W K > 1, we obtain a yield of 16 532 sig- 
nal events, (c) Requiring W n — Wk > 4, we obtain a yield of 10 331 K~K + signal 
events. The vertical dashed lines indicate the signal and sideband regions used for 
the lifetime and yep fits. 

Because the D° — > K~ K + signals have significant reflection backgrounds due 
to misidentified D° — > K~n + decays, we use a modified version of the mass 
sideband subtraction fitting technique used in our preceding experiment [3]. 
We will discuss the complete technique here. We fit the reduced proper time 
histogram for D° — > K~K + or D° — > K~n + signal region events to a cor- 
rected exponential distribution for the signal added to a background reduced 
proper time histogram obtained directly from combinations falling in either a 
high or low mass sideband. Because this technique assumes that the events 
in symmetrically placed mass sidebands have the same lifetime evolution as 
events in the background within the signal region, it must be modified in light 
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of the K~ir + misi identification reflection shown in Figure 1 that only populates 
the upper sideband. In order to accommodate this reflection, we subtract the 
expected contribution from the D° — > K~tt + reflection from the raw upper 
sideband, reduced proper time histogram. The overall normalization of the re- 
flection contribution is computed from the integral over the sideband domain 
of the reflection peak found through the fit to the D° — > K~K + mass spectrum 
illustrated in Figure 1. The reduced proper time histogram shape is computed 
using the lifetime of the D° — > K~ir + . This lifetime is taken from a joint fit 
of the D° K~K + and K~ir + reduced proper time distributions where the 
four fit parameters are (1) lifetime of the D° — > K~tc + , (2) the y C p parameter 
which relates the D° — > K~K + lifetime to the D° — > K~n + lifetime, (3) the 
background level under the K~ir + signal, and (4) the background level under 
the K~ K + signal. The normalization of the exponential time evolution of the 
D° — > K~K + and K~ir + signals is given by the total number of events in the 
signal region minus the background level fit parameters. 

The reduced proper time contributions for the D° — > K"tt + , D° — > K~K + 
and misidentified D° — > K~n + reflection are of the form f(t') exp(— t'/r) 
where fit') represents efficiency and absorption corrections to a pure exponen- 
tial decay with lifetime r. The use of a multiplicative "efficiency" correction, 
rather than an integral over a resolution function is justified since our reduced 
proper time resolution is less than 1/10 th of the D° lifetime. Because of the 
somewhat large (200 fs) bin widths, we actually integrate the exponential over 
the domain of the bin in computing the signal contribution rather than just 
evaluating the exponential at the bin center. 

A separate f(t') correction factor, determined using a Monte Carlo simulation, 
is used for each of the three exponential contributions. Our Monte Carlo sim- 
ulation includes the Pythia model for photon-gluon fusion and incorporates a 
complete simulation at digitization level of all detector and trigger systems, 
with all known multiple scattering and particle absorption effects. We have 
confirmed that it accurately reproduces the momentum and P± distribution 
for D mesons and the multiplicity and profile of the primary vertex. The Monte 
Carlo was run with 20 x the statistics of the experiment. 

Figure 2 shows the efficiency and absorption corrections (f(t')) obtained us- 
ing this Monte Carlo for both decay modes in 200 fs bins of reduced proper 
time. The /(£') function is obtained by dividing the simulated reconstructed 
charm yield in each reduced proper time bin by the input decay exponential 
integrated over the bin. Figure 2 shows that the Monte Carlo corrections are 
typically less than 5% for both decay modes and the corrections for K~K + are 
statistically consistent with those for K~n + . Sources of potential /(£') varia- 
tion include a minor relative depletion at low t' since the charm secondaries 
must lie within the fiducial area of the downstream microstrips, a depletion 
at large t' for charm candidates decaying downstream of the first microstrip 
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detector, and a slight depletion at low t' since the upstream charm daugh- 
ters need to travel through more material before exiting a target segment 
for the ~ 30% of events whose secondary vertex lies in target material. The 
charm daughter absorption effect is partially compensated when one considers 
absorption of the charm particle itself which tends to favor low t' for those 
events produced within the target material. We assume that the charm ab- 
sorption cross section is 1/2 of the cross section for neutrons. Uncertainty in 
the charm cross section should cancel when the two charm decay lifetimes are 
divided to form y CP . 

Because FOCUS uses a segmented target consisting of four 6.75 mm thick 
BeO sections, each separated by 10 mm, many decays occur in air. The charm 
absorption is minimized in this configuration creating only minor corrections 
to the fitted lifetimes. Except for the « 24% difference between the absorption 
cross section for kaons and pions in the momentum range relevant to the D 0, s 
reconstructed in FOCUS, the minute absorption correction will be common 
to both decay modes [4] and cancel in yep ■ To give a feeling for the smallness 
of the absorption correction for lifetimes in FOCUS, we offer the following 
example. Computing the fit') correction using a Monte Carlo where the ab- 
sorption cross sections for charm secondaries have been scaled by a factor of 
60% relative to their known values, causes the lifetimes for the K~ K + and 
K~7i + to decrease by only about 1.2 fs which implies a very small absorption 
systematic uncertainty for yep compared to the size of the statistical error. 

The background levels in the D° — > K~ir + and D° — > K~K + signal region 
are parameters in the lifetime fit. We have employed two ways of determining 
these parameters that are used to normalize the background contribution to 
the reduced proper time histogram in the signal region. The first method 
determines the background levels by finding the number which best fits the 
time evolution in the signal region. The second method combines information 
from the lifetime evolution with additional information from the fits to the 
mass distribution such as those shown in Figure 1. We accomplish this by 
adding additional likelihood terms which tend to "tie" the total background 
level to the one deduced from the mass fit. Specifically for the D° — > K~tt + 
we add the log likelihood of a Poisson distribution which ties the D° — > K~ir + 
level to half of the sum of the number of candidates in the upper and lower 
sideband. For the D° — > K~ K + we add a x 2 -like likelihood penalty term 
which ties the background level to the integral of the polynomial used to 
represent the background under the D° — > K~ K + peak in the signal region. 
The incorporation of information from the mass fit tends to reduce errors by 
15-20% compared to the fits where the background level is determined from 
the time evolution alone. 

Figure 3 shows the t' evolution for the D° — > K~n + and D° — > K~K + along 
with the predicted number from the lifetime fit. The confidence level for these 
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fits are 2% and 55%, for the K tt + and K K + , respectively. From these fits 
we obtain a lifetime asymmetry of: 



T(CP even) + T(CP odd) 



and a D° — > K tt + lifetime of 



71 



409.2 ± 1.3 fs (statistical error only). 



Using our value for the fitted lifetime asymmetry, we compute t(D° — K K + ) 
395.7 ± 5.5 fs (statistical error only). 
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Fig. 2. Monte Carlo correction factors for D° — > K~tt + and K~K + for Ija > 5 
and W n — Wk > 4. We have offset the K~K + points slightly for clarity and have 
given them "flats" on their error bars. Monte Carlo corrections are rather slight 
with these cuts and the corrections for D° — > K~7r + are the same within errors as 
those for D° -> K~K+. 



The systematic uncertainty was studied by gauging the variation of the fitted 
lifetime estimates as the analysis cuts and fitting technique are varied. Because 
the Monte Carlo corrections shown in the f(t') plots (Figure 2) are consistent 
between the D° — > K~7i + and K~K + sample, one expects the dominant 
systematic error on yep to come from potential differences in the background 
under the K~K + peak. By varying the minimum Ija cut from 5 to 9 for the 
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Fig. 3. Signal versus reduced proper time for D — ► K~ir + and K~ K + requiring 
Wk-W k >4 and £/a > 5. The fit is over 20 bins of 200 fs bin width. The data is 
background subtracted and includes the (very small) Monte Carlo correction. 




Fig. 4. Stability of the yep results for 9 sets of clean-up cuts and 4 different fit op- 
tions. The set of 9 cut variants consists of 3 different kaonicity cuts (AWk > 1, 2, 4), 
each with three different detachment cuts if a > 5, 7, 9. The first 18 values use a 15 
bin fit; the last 18 values use a 20 bin fit, where the bin size remains fixed at 200 
fs. The 1st and 3rd set of 9 values obtain the background level entirely through the 
time fit. The 2nd and 4th set use the background level which incorporates additional 
information from the mass fits shown in Figure 1. The RMS spread in these values 
is 0.63 % which is considerably smaller than our statistical errors. The horizontal 
lines give our quoted result on yep plus or minus its quoted statistical error. 

case of K~tt + and K~K + , we significantly change the relative background 
level by eliminating non-charm backgrounds. Denning the signal-to-noise as 
the ratio of the signal height to background height at the location of the 
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Fig. 5. Stability of the r(if _ 7r + ) results for 9 sets of clean-up cuts and 4 different 
fit options. The 36 estimates are plotted according to the convention of Figure 4. 
The RMS spread in these values is 0.28 fs which is considerably smaller than our 
statistical errors. The horizontal lines give our quoted result on t(D° — » K~n + ) 
plus or minus its quoted statistical error. 



Gaussian signal peak, the fits used to measure the lifetimes have signal-to- 
noise ratios which range from 8.9 to 19.3 for the K~tt + , from 2.3 to 5.7 for the 
K~K + . Changing the Cerenkov log likelihood ratio cuts from AWk = 1 to 4 
significantly changes the level of charm reflection contamination (as evidenced 
by Figure 1), and reduces contamination from combinatoric background. 

Figure 4 demonstrates the stability of the yep results for 9 sets of clean- 
up cuts and 4 different fit options. The 9 clean-up cut variants considered 
were 3 different kaonicity cuts (AWk > 1,2,4), each with three different 
detachment cuts Ija > 5, 7, 9. Each set of 9 points in Figure 4 are ordered as 
(AW K > 1, i/<r > 5), (AW K > 2, ijo > 5), (AW K > 4, ija > 5) , (AW K > 1, 
i/a > 7), (AW K > 2, i/a > 7), (AW K > 4, i/a > 7), (AW K > 1, £/a > 9), 
(AW K > 2, ija > 9), and (AW K > 4, ija > 9). The four fit option results, 
summarized in Figure 4, include varying the lifetime range and the method 
used to obtain the background normalization. We show the fitted D° — > K~n + 
lifetime for each of the 9 cut variants and 4 fit options in Figure 5. Our 
quoted systematic error was evaluated by first calculating the shifts in y CP for 
three different detachment cuts, three different kaonicity cuts, two different 
background normalization options, and three different lifetime fit ranges (10, 
15, and 20 bins of 200 fs). These shifts were then combined in a conservative 
manner by adding them in quadrature to obtain the quoted systematic error. 

Additional studies, beyond those used to obtain the quoted systematic er- 
ror, were made to assess the validity of the Monte Carlo corrections and the 
background subtraction technique. We studied the consistency of the K~n + 
lifetime from four samples split according to momentum and primary ver- 
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tex location for a variety of detachment and Cerenkov cuts. These studies 
include comparing the K~n + lifetime for 4 samples split according to the 
K~n + momentum and the location of the primary vertex (in the upstream 
two targets versus downstream or upstream two targets). Although the Monte 
Carlo (fit')) corrections for the high momentum, downstream sample were 
much severe than those for the other three samples , the K~n + lifetimes for 
all 4 samples were consistent within errors. 

To assess the validity of the background subtraction method, we fit the ycv 
values using two versions of the same set of 36 fit variants summarized in 
Figure 4. The first version used sideband regions that were set to 1/2 of the 
width of those shown in Figure 1. The second version restricted events to the 
tagged path only (but with full sideband width illustrated by Figure 1). The 
signal to background in tagged-only sample was 4 times larger than that for our 
standard tagged and inclusive combined sample. Besides reducing the overall 
background level, the tagging requirement should have substantially reduced 
any potential backgrounds from partially reconstructed and Cerenkov mis- 
identified D + or A+ decays that have a substantially different lifetime than 
the D° . Significant variations in the result when the sideband width is changed 
might signify the presence of possible backgrounds (such as those from possible 
charm reflections) which have a mass spectrum which is sufficiently non-linear 
that the average of the high and low sideband reduced proper time evolution is 
no longer an inadequate representation for the t' evolution of the background 
in the signal region. 

No problems with the background technique were uncovered by either of these 
alternative versions of the 36 standard fits. In all cases the yep value was 
consistent with our standard value. 80% of the 36 "tagged only" fits and 94% 
of the 1/2 sideband fits returned yep values within the dashed horizontal lines 
of Figure 4. The remaining fits were at most 3/4 of their error bar away from 
these lines. 

We have presented a new measurement of the lifetime ratio between a CP 
even final state, D° — > K~K + and a CP mixed decay, D° — > K~n + of yep = 
(3.42 ±1.39 ±0.74)%. Our analysis techniques have been designed to minimize 
the relative systematic errors between these samples, rather than to obtain 
the best statistical error on the D° lifetime under the assumption of a pure 
exponential decay. For example, inclusion of the D° — > K~7t + 7t~tt + decay 
mode would essentially double our statistics for the D° lifetime. In addition, 
there are systematic error sources such as the overall distance scale error which 
affects our absolute lifetime but not the lifetime ratio which is the principle 
result reported here. 

E791 [1] measures Ar = 2(r KK -T Kn ) = 0.04 ± 0.14 ± 0.05 ps _1 . Combining 
this and their measurements of the KK and Kn lifetimes, we obtain a value 
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of yep — (0-8 ± 2.9 ± 1.0)% which is consistent with our measured value of 
y CP = (3.42 ± 1.39 ±0.74)%. 

A more recent result exists from the CLEO Collaboration. CLEO searches 
for mixing effects by studying the possible interference of mixing with direct 
doubly-Cabibbo-suppressed decays in the time evolution of D* + — > 7t + {K + tt~) 
decays [5]. They report a 95% confidence level range on a variable they call y' of 
—5.8% < y' < 1%. If the level of CP violation in charm decays is negligible, the 
CLEO y' variable is a rotational transformation of the yep variable reported 
here and a variable which depends on the CP eigenstate mass difference with 
the angle of rotation being due to a strong phase shift. Theoretical estimates 
on the size of this angle differ significantly [6] making a precise comparison of 
our result with the CLEO result impossible at the present time. 

Because of our high statistics, the error on yep reported here can be reliably 
interpreted as a Gaussian error for the purposes of combining with other mea- 
surements {e.g. The \ 2 versus y^p in a ±1 a domain about the fit minimum 
is well fit by a parabola). Our measurement represents the most precise di- 
rect measurement of the neutral D meson CP eigenstate lifetime difference. A 
wide range of Standard Model and non-Standard Model predictions on mixing 
through the width difference appear in the literature [7]. 

We wish to acknowledge the assistance of the staffs of Fermi National Ac- 
celerator Laboratory, the INFN of Italy, and the physics departments of the 
collaborating institutions. This research was supported in part by the U. S. 
National Science Foundation, the U. S. Department of Energy, the Italian Isti- 
tuto Nazionale di Fisica Nucleare and Ministero dell'Universita e della Ricerca 
Scientifica e Tecnologica, the Brazilian Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico, CONACyT-Mexico, the Korean Ministry of Educa- 
tion, and the Korean Science and Engineering Foundation. 



14 



References 



[1] E791 Collab., E. M. Aitala et al, Phys. Rev. Lett. 83 (1999) 32. 

[2] E687 Collab., P. L. Frabetti et al, Nucl. Instrum. Methods Phys. A 320 (1992) 
519. 

[3] See for example E687 Collab., P. L. Frabetti et al Phys. Rev. Lett. 70 (1993) 
1755 ; Phys. Rev. Lett. 71 (1993) 827 ; Phys. Lett. B323 (1994) 459. 

[4] Particle Data Group, C. Caso et al, Eur. Phys. J. C3 (1998) 1. 

[5] CLEO Collab., R. Godang et al, CLNS-99- 1659, Dec 1999, submitted to Phys. 
Rev. Lett. e-Print Archive: |hep-ex/OOOT060| . 



[6] A. F. Falk, Y. Mr, A. A. Petrov, JHEP 9912:019,1999. 



[7] Harry Nelson, UCSB HEP 99-08 (Aug. 1999), e-Print Archive: |hep-ex/00010p0| 



and references therein; E. Golowich, Proceedings of the Conference on B Physics 
and CP Violation, Honolulu, HI (March 1997) and references therein; A. A. 
Petrov, Phys. Rev. D56 (1997)1685; E. Golowich and A. A. Petrov, Phys. Lett. 
B427 (1998) 172; J.F. Donoghue, E. Golowich,B.R. Holstein and J. Trampetic, 
Phys. Rev. D33 (1986) 179; T. Ohl et al, Nucl. Phys. B403 (1993) 605; S. 
Bergmann, Y. Grossman, Z. Ligeti, Y. Nir, A. A. Petrov , FERMILAB-Pub- 
00/102 (May 2000), e-Print Archive: [hep-ph/00051~81 . 



15 



